NASA Technical Memorandum 87004 


Vibration and Buckling of Rotating 
Pretwisted, Preconed Beams 
Including Coriolis Effects 


(liSA-Tfl-67004) VIEBATION ANE BUCKLING OF 
B01AIING, FBBIbISlIE, FEBCCSIE BliHS 
IBCLUCING C00EICLIS 1IEEC1S (NASA) 21 p 
HC A3 2/HF A01 CSCI 20K 

G3/39 


K.B. Subrahmanyam and K.R.V. Kaza 
Lewis Research Center 
Cleveland , Ohio 


Prepared for the 

10th Biennial Design Engineering Conference and Exhibit 
on Mechanical Vibration and Noise 
sponsored by the American Society of Mechanical Engineers 
Cincinnati, Ohio, September 10-13, 1985 


NASA 


N85-258S3 

Unclas 

21133 


E- 231 0 


VIBRATION AND BUCKLING OF ROTATING, PRETWISTED, PRECONED BEAMS INCLUDING CORIOLIS EFFECTS 


K.6. Subrahmanyam 1 and K.R.V. Kaza 
National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135 


ABSTRACT 

The effects of pretwist, precone, setting angle 
ana Coriolis forces on the vibration and buckling 
behavior of rotating, torsional ly rigid, cantilevered 
beams are studied in this investigation. The beam is 
considered to be clamped on the axis of rotation in 
one case, and off the axis of rotation in the other. 
Two methods are employed for the solution of the 
vibration problem: one based upon a finite-difference 
approach using second order central differences for 
solution of the equations of motion, ano the other 
based upon the minimum of the total potential energy 
functional with a kitz type of solution procedure 
making use of complex forms of shape functions for the 
oepenaent variables. Numerical results obtained by 
using these methods are compared to those existing in 
the literature for specialized simple cases. Results 
indicating the individual and collective effects of 
pretwist, precone, setting angle, thickness ratio and 
Coriolis forces on the natural frequencies and the 
buckling boundaries are presented and discussed. 
Furthermore, it is shown that the inclusion of 
Coriolis effects is necessary for blades of moderate 
to large thickness ratios while these effects are not 
so important for small thickness ratio blades. 

Finally, the results show the possibility of buckling 
due to centrifugal softening terms for large values of 
precone and rotation, 

INTRODUCTION 

An important phase in the development of advanced 
turboprop blades currently in progress at the Lewis 
Research Center Is the development of analytical blade 
models that are capable of predicting the vibration 
ana flutter characteristics to an acceptable degree of 
accuracy. The turboprop blades are of thin cross 
sections with large, variable sweep and pretwist, and 
are mounted on a rotating hub at a setting angle. 
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Moreover, the blades are subjected to considerable 
centrifugal loading and develop steady-state 
deflections that are large compared to the 
cross-sectional dimensions. These blades, under 
certain conditions are subjected to centrifugal 
compressive forces on certain cross sections, and 
consequently the possibility of buckling cannot be 
ruled out. The recent research in the vibration and 
flutter of these blades indicate that the effects of 
nonlinearities, variable sweep, and pretwist must be 
included in the analysis for accurate prediction of 
the blade dynamic characteristics. The finite-element 
modelling of these blades is appropriate for the 
determination of their dynamic characteristics. 
However, such studies with the existing codes showed 
that the predicted results are satisfactory only for 
the first few mooes and that the complicating effects 
included in these theories make the understanding of 
the individual and collective effects of the governing 
parameters impossible. These considerations have lead 
to the modelling of the turboprop blade configurations 
using the simpler beam theories with the complicating 
effects successively taken into account to establish a 
physical understanding and to reveal the relative 
importance of the individual effects. A preliminary 
part of such : ly using a simplified beam model with 
a set of lin?a vquations of motion is reported in 
this paper. In the beam model the effects of sweep on 
the dynamic behavior of the blade are introduced by 
setting the axis of the blade at an angle with respect 
to the plane of rotation. This angle is usually 
called the precone angle. The more complicated 
nonlinear aspects form the subject matter of a 
subsequent paper. 

An examination of the existing literature using 
beam theories indicates that several aspects of the 
vibration and buckling of rotating beams mounted off 
the axis of rotation were considered by several 
investigators (W)* The effects of Coriolis forces 
on the vibration of radial rotating beams were 
considered in (8,9) among others. In (5), an axially 
prestresseo, pretwisted rotating blade mounted off the 
axis of rotation was considered and an extended 
Galerkin procedure was used for solution of the 
equations of motion. In all other Investigations 


cited, the beams were considered to be untwisted. 
Reference (J.) uses a method of successive 
approximations, (2,3) use a perturbation technique 
while (2) gives asymptotic expansion formulas together 
with a composite expansion formula for determination 
of the buckling boundary of untwisted rotating blades 
mounted off the axis of rotation. Reference (4) uses 
an integrating matrix approach for the solution of 
untwisted blade stability analysis, while (6) unifies 
the rotating beam vibration analysis allowing for 
precone and second degree non linearities. From the 
published literature, it appears that the combined 
effects of pretwist, large precone, setting angle and 
Coriolis effects on the vibration and buckling 
behavior were not stuoied. Furthermore, greater 
emphasis had been laid on the solution of the 
governing equations of motion using various methods of 
solution while solutions using energy functionals have 
received considerably less attention. 

The specific objectives of the present paper 
are: (1) to solve the pretwisted, preconed rotating 
cantilever blaae vibration problem allowing for the 
Coriolis effects for two cases of root clamping, that 
Is, on the axis, and off the axis of rotation, and (2) 
to study the effects of pretwist, precone, setting 
angle, Coriolis forces, and the ratio of the principal 
moments of inertia on the vibration frequencies and 
buckling boundaries. Two methods of solution are 
employed: (a) a finite difference procedure using 

second order central differences for the solution of 
the governing equations of motion and, (b) the total 
potential energy functional ana the Rltz type solution 
using complex shape functions for the dependent 
variables. The present solution methods provide a 
positive check on the accuracy of the results for 
untwisted blade cases in that the finite difference 
procedure gives close lower bound solutions while the 
potential energy method gives upper bound results. 
Further validation of the results is accomplished by 
comparing the present results to closed form exact 
solutions wherever possible. The solution methods and 
parametric results provide a sufficient insight for 
further extensions of the present effort to account 
for blade sweep and second degree nonlinearities in 
the analysis. 

FREQUENCY EQUATION BY SECOND ORDER FINITE-DIFFERENCE 
METHOD 

The linear set of coupled bending-bending 
equations of motion for a torslonally rigid slender 
beam (Fig, 1) of rectangular cross section allowing 
for the effects of linear pretwist, precone and 
Coriolis forces can be shown (10) to be (a list of 
notation is given in Appendix EJ: 
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It may be noted that the last two terms in 


Eps. (1) and (2) are shown with */* signs. For a 
cantilever blade case with root clamping on the 
periphery of a disc of radius R ana the blade 
protruoing radially outwards, the upper sign is to be 
taken. For the case with root clamping on the rim of 
a ring of radius R and the blade protruding radially 
inwards, the lower sign Is to be assumed. 

It is of Interest to examine the coupled 
bending-bending equations and point out the various 
important terms arising due to the Inclusion of 
precone. It may be seen that in Eq. (1), the term 
"Pc) appears due to the inclusion of 
preconing which is a softening term. Even for the 
case of a rotating beam protruding radially outwards 
and under centrifugal tension, the effect of this term 
is to reduce the effective stiffness, thereby 
provlolng the mechanism for Instability at reasonably 
large values of precone and rotational speed. It may 
also be noted that the Coriolis terms, which are 
underlined In Eqs. (1) ana (2), are linear In v and 
w and are associated with precone. In the absence of 
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precone, the equations do not contain linear Coriolis 
forces. As could be seen in the next sections, the 
effect of these Coriolis terms is extremely important 
for blades of moderate to large thickness ratios, and 
also for blades possessing large pretwists. 

In seeking a solution by using a 
finite-difference procedure, one substitutes the 
finite-difference expressions for the derivatives into 
the Differential equations, for any arbitrary station 
j inside the beam domain. When the beam is divided 
into n segments, j varies from 0 , 1 , 2 ,.. to n; the 
clamped end being denoted by 0 and the free end by n. 
A set of n equations are written by successively 
assigning values for j from 1 to n. As can be seen 
from the expressions for the derivatives of an 
arbitrary variable t at the station j in terms 
of the second order central differences given in 
Appendix A, one has to eliminate the evaluation of the 
functions at the fictitious stations outside the beam 

namely w_ 2 , w_i, *n+l» w n+2* w c + 3* v -2* 
vi, v Q +i, Vfl+ 2 * ana Vp+ 3 . This can be 
accomplished by using the boundary conditions together 
with the recursive relations developed and discussed 
in detail in (_U,12). These relations are also 
included in Appendix A for ready reference. One can 
thus write a set of 2 n equations in terms of the 

variables W|, w?...w n , v l* v 2 «** v n 

the coupled bending-bending equations of motion given 

by Eqs. ( 1 ) ana ( 2 ) for values of j - l, 2 ,...n. 

These equations can be represented in the following 
matrix form: 
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where ml, m2. Cl, kl, k2, k3 and kj are square 
matricerof "order iT x n~ each, {w} and 
(v) are column matrices containing the vector 
wi, W 2 .*.w n ; vi, V 9 *«'V n and £ is null 
matrix. For brevity elements of these matrices are 
not presented here. Equation (4) is of the form 
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The eigenvalues of Eq. (9) for a conservative system 
are real and exist in pairs Pi, -pi, P?* 

-P2,...p n , -p n ; n being the order of M, £ 
and K. However, eigenvectors for Eqs. ( 6 ) and (9) 
are the same and will occur as pairs of complex 
conjugate vectors (13) . For instabilities, the 
frequency p reduces to zero (static buckling). 

FREQUENCY EQUATION BY USING ENERGY FUNCTIONALS 

The Lagrangian functional for a rotating, 
pretwisted, preconed blade including Coriolis effects 
but disregrading nonlinear effects can be derived 
from ( 6 ,_ 1 £) and is given by 
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In Eq. (11), the last term is shown with T /± 
signs as before with the convention that the upper 
sign is to be taken for the cantilever blade case 
with root clamping on the periphery of a disc and the 
blade protruding outwards, while the lower sign is to 
be assumed for the case with root clamping on the rim 
of a ring and the blade protruding radially inwards. 

In order to formulate the frequency equation 
through the Ritz process, shape functions for the 
variables w and v are assumed in the following 
complex form since the out-of -phase Coriolis effects 
are also included in the analysis (9): 
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j J 
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where A^ , Bi, Cj, Dj are the undetermined 
arbitrary real parameters in the shape functions, 
fj is the polynomial function defined as 
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The shape functions assumed here satisfy the boundary 
conditions 

v«w*v'=w' *0 at n * 0 and 

v' , «w , '=v ,, ‘*w ,,, *0 at n«l (15) 

The real and imaginary parts of the shape functions 
are substituted into the Lagrangian functional 
separately, the indicated integrations are performed, 
the resulting functional is time averaged in the 
traditional manner and the Ritz process i^ applied by 
minimizing the resulting Lagrangian functional, that 
is 

aL aL aL aL 
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The resulting equations can be written in the 
following form of the eigenvalue problem after 
transformations similar to those oescribed in the 
previous section are maoe: 

C + pD - 0 (17) 

Where £ ana £ are complex matrices. 

RESULTS AND DISCUSSION 

The eigenvalue problems defined by Eqs. (9) and 
(17) were solved by using computer programs developed 
in FORTRAN language. These programs were run on IBM 
and CRAY computers at the NASA Lewis Research 
Center. Use was made of standard eigenvalue 
extraction routines (14), ana a 15-point Gaussian 
quadrature formula for carrying out integrations 
where necessary. Results thus generated are 
presented and discussed in what follows. 

CONVERGENCE 

The relative convergence rates produced by the 
finite-difference method using second order central 
differences, and the potential energy method were 
obtained for several configurations of rotating blade 
cases. A typical set of such results are presented 
in Table I for the case of a preconed, rotating beam, 
including Coriolis effects. From the convergence 
pattern shown in Table I, it can be seen that the 
potential energy method produces convergence to seven 
significant figures with an eight-term solution 
(k - 8) while the finite-difference method produces 
the lowest five frequencies to within 0.05 percent 
error when the beam is divided into 30 segments 
(n ■ 30), although the convergence is up to four 
significant figures. It is worth noting here that 
the convergence resulting from the present 
formulation of the finite-difference procedure making 
use of second order central differences is much more 
rapid than a solution that could be obtained by the 
classical finite-difference method. It has been 
established Q2,j5) that this second order 
finite-difference method is computationally more 
efficient than the classical approaches of solving 
the equations of motion using ordinary 
finite-difference procedures. Further, the 
Richardson's extrapolation procedure (16,17) can be 
applied for frequencies obtained with Two or three 
different matrix sizes to Improve the accuracy while 
keeping the matrix sizes small In the 
finite-difference methods. 


When a comparison of the present converged 
results is made to those presented by Leissa and Co 
(JO for the specialized simple case presented in 
Table I, it is observed that present results show 
excellent agreement with those of (9). The Ritz 
method of solution was employed in T9) making use of 
a different form of shape functions. Further, the 
extensional aeformantion was taken into account in 
(_9) while the present work treats the blades to be 
extensional ly rigid. It is thus evident from the 
close agreement of the present results with those of 
(9) that for the geometric and physical parameters 
considered in this work, the extensional deformation 
can safely be ignored. Exclusion of this additional 
degree of freedom reduces the number of equations of 
motion ana the consequent computational space and 
time. 

Finally, for all the vibration and stability 
problems reported in this work, solutions were 
obtained by using the upper bound potential energy 
method with k * 8, and the usually lower bound 
finite-difference method with n - 30. In view of 
the close agreement between these two sets of results 
observed (of the order of 0.05 percent difference), 
only one set of results are presented in what follows 
for brevity. 

VIBRATION AND STABILITY OF ROTATING BLADES 

In order to study the effects of pretwist, 
precone, Coriolis forces and thickness ratio (d/b) on 
the natural frequencies of rotating blades, several 
cases of twisted, rotating blades were solved. 

Typical cases of thickness ratio that approximately 
represent propeller blades (d/b = 0.5 or 0.25), and 
advanced turboprop blades (d/b «= 0.05) were 
considered for a typical aspect ratio (L/b) of 5. 

The dimensionless rotational speed (q/u>i) 

(where a and wi are respectively the 
rotational speed, and the uncoupled fundamental 
bending frequency of the blade) was varied. Typical 
set of such results are presented in Table II to IV, 
and in Figs. 2 and 3. Table V is prepared using the 
results presented in Tables III and IV to Illustrate 
the effect of Coriolis forces on a typical preconed 
blaae. It may be noted that K is set equal to zero 
in these calculations for simplicity. 

The effect of pretwist on nonrotating blade 
frequencies predicted by beam theories is well 
understood in the published literature. The 
frequencies presented in Table II for ft/u»i ■ 0 
case agree well with those existing in the 
literature, and exhibit the well established coupling 
trends (18-20). The fundamental mode frequency is 
seen to Increase slightly with increasing pretwist, 
consistent with the Inherent quality of the beam 
theory (J21). However, for rotating pretwisted blades 
it can be seen that the fundamental mode frequency 
shows a decreasing trend for increasing pretwist 
beyond certain values of a /»!• 

The effect of rotation on the coupling trends of 
pretwisted blades is illustrated in Table II. For 
pretwisted blades with thickness ratio greater than 
0.1596, the second coupled mode will be closer to the 
fundamental uncoupled mode in chordwise direction 
while the third coupled mode will be closer to the 
uncoupled second mode frequency in the flapwlse 
direction. The effect of pretwist for such blades Is 
to couple the second mode uncoupled flanwlse 
frequency ana the fundamental chordwise frequency in 
such a way that the lower frequency is reduced and 
the higher one is Increased with Increasing 
pretwists. This important coupling t end observed 
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for nonrotating pretwisted blades does not 
necessarily bold for rotating blades. This may be 
verified from the results presented in Table II for 
pretwisted blaoe cases with thickness ratio of 0.5 
and rotational speed parameter values of 3 and 
greater. The reason for the change in the coupling 
trend is th^t although the coupled frequencies are 
raised by the rotation induced centrifugal forces, 
not all the modes are stiffened in a proportional 
manner. 

Next, the effects of pretwist and precone on 
rotating blaae frequencies are considered. Some 
typical results are presented in Table III neglecting 
Coriolis effects, and in Table IV including Coriolis 
effects. Futher results presented in Table V 
indicate the effect of Coriolis forces, those in 
Fig. 2 indicate the effect of pretwist, and those in 
Fig, 3 indicate the effect of precone in a parametric 
form. The following observations are made from the 
results presented: 

(1) For small precone angles and small 
pretwists, an increase in the rotational speed 
increases the frequencies of all modes. For large 
precone angles, the frequencies of pretwisted or 
untwisted blade Increase in a nonlinear manner up to 
a certain value of the rotational parameter, 
q/w 1 , beyond which a decreasing trend is 
observed. At some value of the rotational parameter 
a/wl, for a given thickness ratio and precone, 
the fundamental mode frequency becomes zero 
indicating static instability. A aetailed discussion 
on this phenomenon for untwisted blades was presented 
in (6,22). This instability occurs over a zone of 
n/a>i values indicating an instability bano, 
beyond which the blade becomes stable again. This 
can be verified by referring to Table IV for precone 
angles of 60 and 90 ov large thickness ratio blade 
cases. By comparing the results with and without 
Coriolis effects presented in Tables III and IV 
respectively, it can be seen that the Coriolis 
effects change the onset and width of the instability 
band for large thickness ratio blade cases. However, 
the Coriolis effects appear to be not so significant 
for small thickness ratio blade cases operating at 
low rotational speeds. 

(2) The effect of precone in the absence of 
Coriolis effects can be seen by comparing the results 
presented in Tables II and III. From these tables, 
it can be seen that the frequencies of all modes 
decrease as the precone is increased, at a given 
rotational speed for all thickness ratios and 
pretwists. For smaller rotational speeds, the blades 
are stable for all precone angles. For n/«i 
values of one and greater, instabilities could occur 
over a zone of precone angles depending on the 
pretwist and thickness ratio. However, inclusion of 
Coriolis effects makes these frequency va» atlons 
somewhat different. This fact can be verified by a 
mutual comparison of the results presented in 
Tables III and IV for a 30* pretwisted blade with 
fl/«l « 3. It can be seen that for large 
thickness ratio blades with d/b . 0.5, the 
instability band starts for precone angles greater 
than 45 , and extends up to 90* when the Coriolis 
effects are Ignored while the instability band Is 
restricted to between 47* and 50* when the Coriolis 
effects are included. 

(3) In order to provide a clearer Insight into 
the Coriolis effects, the percentage variations of 
the frequency ratios were calculated from the results 
presented In Tables III and IV for a typical precone 
angle of 45 . These results are presented in 
Table V. It can be seen from these results that the 


effect of Coriolis forces becomes more significant as 
the rotational speed is increased. The fundamental 
mode frequency shows a decreasing trena when the 
Coriolis forces are included while the higher mode 
frequencies show both increasing and decreasing 
trends. For large thickness ratio blades, the 
Coriolis effects produce large frequency ratio 
variations even for low rotational speeds. 

Significant frequency ratio variations occur for the 
first and the second inodes of untwisted blades, and 
even for higher modes for pretwisted blades. The 
Coriolis effects are important even for low thickness 
ratio blades for high rotational speeds, particularly 
when pretwist is present. 

(d) The effect of varying pretwist for a given 
thickness ratio and precone angle are presented in 
rig. 2. Specifically, for ap c « 45* and 
thickness ratio of 0.25, the semi-indefinite 
instability band starts at a o/«i value 
somewhere between 5.7 to 5.8 for a 60 pretwisted 
blade while such an instability band starts at a 
fi/wi value somewhere between 5.1 to 5.2 for a 
SO pretwisted blade. For small thickness ratio 
blades with d/b ■ 0.05 and ap c * 45*, instability 
starts at 27.5 < a/un < 28.0 for y - 60*, and 
at 23.0 < a/u>i < 23.5 For y * 90*. For 
fPc . 60 , d/b * 0.05, instability starts at 
1.48 < uUi < 1.49 for blades of 0* or 30* 
pretwTst. For e Pq « 90*, d/b - 0.05 and y - 0* 
blade case, instability starts at n/ w i « 1.0, 
while for the same blade with 30* pretwist, 
instability occurs at 1.003 <a/u»i < 1.004. 

These observations reveal that an increase in 
pretwist can either increase or decrease the 
rotational speeo to cause instability depending on 
whether the precone is lesser or greater than 45* 
respectively. . 

(5) The effect of Increasing the precone angle 
for low thickness ratio pretwisted blade case is 
shown in Fig. 3. It can be seen from this figure 
that considerable change in the slopes of the curves 
occurs for precone angles greeter than 45*. 

From the foregoing discussion, it may be 
concluded that the Coriolis effects associated with 
precone angle may be neglected for low thickness 
ratio blades subjected to low rotational speeds. 

Since the aovanced turboprop blades possess low 
thickness ratios, ana operate around speeds of the 
order of q lv\ ■ 1, it appears logical to 
disregard the Coriolis effects associated with 
precone. However, it Is essential to retain the 
Coriolis effects when pretwists of larger magnitude 
are present for the blades operating even at low 
rotational speeds. 

BUCKLING OF ROTATING BLADES MOUNTED OFF THE AXIS OF 
R0TA1 ION 


A geometric arrangement giving rise to 
rotatlonally induced radial forces which are 
compressive rather than tensile Is shown in Fig. 4. 
Such rotating beams without pretwist and precone have 
been analyzed for g specialized cases of Inplane 
vibration (9 ■ 90 ) or out-of-plane vibrations 
(? - 0 ) extensively Q-7). For a beam even 
without pretwist and precone, the flapwise and 
lagwise bending motions are coupled for setting 
angles other than 0* and 90*. The purpose of this 
section is to re-examine the buckling behavior of 
rotating beams In the presence of pretwist, precone 
and Coriolis forces using the coupled bending-bending 
equations given earlier. 

An examination of the published work reveals 
that most of the problems in the different analytical 


5 


approaches were due to some convergence problems of 
the solution procedure adopted. These were overcome 
by alternative solution methods (4) or by mating use 
of certain expansion formulae (7J* To ascertain that 
the present numerical procedures do not exhibit 
similar drawbacks as were observed in the existing 
works, specialized simple cases were solved first, 
comparisons of present results were made to existing 
stanaaro values, and accuracies of present 
theoretical developments were established. As a 
typical example, the closed form solution obtained in 
( 2 ) for the case of buckling out of the plane of 
rotation (9 - y » ep c ■ 0 ) gives a value of 

* 5. 6746 for R ». i.O. Present solution with the 
potential energy method gives the values of n/*i for 
R « 1 as 5.6912, 5.6753, ana 5.6747 for k « 2 , 3 
and 4, respectively. For k * 8 or 9, this value 
converges to 5.67468. The finite-difference 
procedure produces a value which is about 0.06 
percent off this exact value. Further, the buckling 
boundaries for in-plane and out-of-plane vibration 
cases plotted in Fig. 5(a) agree very closely with 
those obtained by using the integrating matrix 
methoo ^4). 

A parametric study was made to study the effects 
of pretwist, setting angle and thickness ratio (d/b) 
on the buckling boundary. Results of this 
investigation are presented in graphical form in 
Figs. 5 (b) to (d) for d/b » 0.25, and in Figs. 6 (a) 
to (c) for d/b - 0.05. From the results presented, 
the following observations can be made: 

1. When the blade thickness ratio is small, or 
the ratio of the principal second moments of area is 
large, the^stabil it^ boundaries of pretwisted blades 
for 9 * 0 * anci 90 follow closely those for 
untwisted blaaes. 

2. An increase in the pretwist angle for a given 
value of setting angle 9 moves the stability 
boundary towards left. That is, for a 

given rotational speed, buckling occurs at a lower 
value of R for increasing t ^t a given 9. 

3. An increase in the setting angle 9 For a 
given pretwist has similar effect as discussed in (2) 
above. 

4. For nonzero pretwist, all curves tend to meet 
at R ■ 0 as the rotational parameter, n/Xi, 

tends to infinity. 

5. For large values of d/b, there is a marked # 
change in the stability boundary. ^Comparing 9 « 0 , 
y 15* curve with 9 * 0*, y - 0* curve presented 
in Fig. 5(a), one observes that the former closely 
follows t.he latter up to a low value 

n/x x and detaches from this path to rapidly 
reduce the value of R for increasing values of 
As the pretwist angle is increased, this phenomenon 
becomes more pronounced. 

The effects of precone, pretwist and Coriolis 
forces on the critical speed for causing instability 
are presented in Tables VI and VII. From these 
limited sets of results, it appears that the critical 
speed for causing buckling for a given value R 
decreases very rapidly as the precone is increased 
for values of R up to about 1.2. Beyond this value 
of R, there is a slight gain in the critical speed to 
cause instability as the precone is Increased from 
0*. The Inclusion of Coriolis effects in the analysis 
decreases the value of r at a given rotational 
speed to cause Instability generally. As this 
critical speed is Increased further, the effect of 
the Coriolis terms decreases. The nonlinear effects 
may be more Important in such cases rather than the 
Coriolis terms. The effect of nonlinear terms were 
studied in (6) for simple cases of preconed beam 


vibration. It was observed that approximately 
6 percent decrease in the frequency value occurs for 
first flapwise benaing frequency, and about 
30 percent increase in the first lagwise frequency 
for an untwisted, 15* preconed beam with d/b - 1 and 
n/a>X ■ 3.2 due to the steady state deformation 
induced by the nonlinear terms. From these 
observations, it appears that inclusion of nonlinear 
terms and the Coriolis effects into the analysis is 
imperative. Finally, the effect of increasing 
pretwist, shown in Table VII, reveals that 
for buckling to occur at a given critical speed, the 
value of R may either increase or decrease. 

CONCLUSIONS 

Parametric studies were conducted to ascertain 
the individual and collective effects of pretwist, 
precone, setting angle, Coriolis forces and blade 
thickness ratio on the vibration frequencies and 
buckling boundaries of rotating beams. Two methods 
of solution for studying blade vibration and 
stability were used, namely, a finite difference 
procedure based upon second-order central differences 
that usually produces close lower bound solutions, 
and the potential energy method that produces close 
upper bound solutions. Results obtained by using 
these two methods were found to be in excellent 
agreement. Further validation of the present results 
was accomplished by comparisons to results in the 
literature for specialized simple cases. 

The parametric results show that the Inclusion 
of Coriolis effects is necessary for blades of 
moderate to large thickness ratios while these 
effects are not so important for small thickness 
ratio blaaes. Thus, the linear Coriolis terms 
associated with precone may be neglected in the 
dynamic analysis of advanced turboprop blades. 

However, the effect of Coriolis terms in the presence 
of second degree nonlinear terms is yet to be 
assessed. For a given thickness ratio and pretwist, 
an increase in rotational speed has a destabilizing 
effect for large precone angles. For a given 
thickness ratio and aspect ratio, an increase in 
pretwist angle has # a stabilizing effect for precone 
angle less than 45*, and a destabilizing effect for 
precone anqles greater than 45*. For a beam mounted 
off the axis of rotation, the pretwist and setting 
angle have a significant influence on 
rotation-induced buckling instability. However, this 
influence depends markedly on the blade thickness 
ratio. 

The parametric results presented in this work 
are believed to be useful for future comparisons of 
theoretical solutions Including sweep and nonlinear 
effects, and thereby ascertaining the importance of 
these complicating effects. 

APPENDIX A 

FINITE-DIFFERENCE EQUATIONS FOR DERIVATIVES 

The finite-difference expressions for the 
derivatives of an arbitrary functional 9(n) in 
terms of second order central differences are given 
as follows for an arbitrary station j: 

*j " [*j-2 ' ®*j-l + 8 Vl - 12h < A1 > 


6 


* l6 *j-l - 30 *j 


length of each elemental beam segment 


+ 16* 


j + l 


- 


12h 


(A2) 




area moments of inertia about major 
and minor principal centroidal axes, 
respectively 


*j " [*j-3 " *j-2 + 13 *j-l ‘ 13 *j + l 

+ 8 V2“ Va]/ 8 ” 3 < A3 > 

♦j V "foV 12 *j-2 - 39 *j-l + 56 *j 

- 3 Vi + 12 *j+2 - Va]A 4 < A4 > 

Equations (Al) to (A4) are applicable for a beam 
di video into n number of segments of equal length 
and n is the nonoimensional length coordinate 
varying from n ■= 0 at root section to n - 1 at 
free end for a cantilever beam. Thus h - 1/n. 

The boundary conditions for a cantilever beam 
under coupled bending-bending motion are 

v«w = v'=w‘ = 0 at n-0 (A5) 

v * 1 * w 1 1 ■ v ' * 1 * w 1 1 1 ■ 0 at n - 1 (A6) 


J 

k 

L 

L 


n 

P 

R 

R 

t 

v 


w 


dummy index 

number of terms in assumed solution 
length of beam 

time averaged value of the Lagrangian 
functional 

number of beam segments 
natural radian frequency 

radius of disc 
nondimensional radius, R/L 
time 

displacement in the chorowise 
direction 

displacement in the flapwise 
direction 


Using the expressions for the derivatives given 
by Eqs. (Al) to (A3) and invoking the conditions of 
symmetry as discussed in Q5), one can write the 
following recursive relations which satisfy the 
boundary conditions for suitably large n: 

Vl - Vl " 2(w n " Vl ); 

Vl - Vl ‘ 2 < v n - Vl> < A7 > 

w n+2 ” w n-2 " 4 * w n " w n-l^ ; 

V2 - v n-2 ’ 4 < v n “ Vl) < A8 > 

W n + 3 " w n-3 “ 6 < w n " Vl>« 

V n+3 - v n-3 " 6 < v n " Vl) < AS > 

W -1 “ V w -2 " w 2* v _i - v i» v _2 " v 2 (A10) 

Equations (A7) to (A10) can be used to eliminate the 
fictitious stations outside the beam domain. 

APPENDIX B 
NOMENCLATURE 

A area at any section 

£*£»£•&&• etc * matrices 

A j* B j» c j»°j arbitrary parameters In shape 

functions 


x 


p Pc 


Y 

p 

“1 


0 


9 


T) 


*1 

( * ) 

( )' 


running coordinate along longitudinal 
x-axis 

precone angle (The angle the blade 
centerline makes with the plane of 
rotation, see Fig. 1). 

pretwist of the blade over length L 

mass density of blade material 

exact fundamental mode frequency of 
straight beam, 3*51602 xi 

rotational speed, radians/sec 

geometric pitch angle, (* + Y n) 

stagger angle (collective pitch, the 
rotation of the blade root chord 
about the blade axis) 

nondimensional length coordinate, x/L 

nondimensional rotational parameter, 
EI nn /()AL4 ‘ l2 

frequency parameter, yEI nn /pAL 4 

dot over a parameter represents time 
derivative, a/at 

prime denotes differentiation with 
respect to n , a/an 


b,o 

d/b 

E 

f j 


breadth and thickness of blade 
thickness ratio 
Young's modulus 
shape function 
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TABLE I. - CONVERGENCE PATTERN OF FREQUENCY RATIOS (p/Xj) OF A PRECONED, ROTATING BEAM 

INCLUDING CORIOLIS EFFECTS 




[ 

— - 0.5, 
“1 

&pc * 30 , 

d/b = 0.5, 

L/d . 20, 

y s R = 0. 





(a) 

Second order f inite-aifference method (lower bound) 



Setting 

angle, 

9 

Mode 

number 

n « 6 

n . 10 

n - 12 

n * 15 

n b 18 

n - 20 

n - 25 

n « 30 


i 

3.62676 

3.63274 

3.63349 

3.63400 

3.63423 

3.63435 

3.63446 

3.63453 


2 

7.29213 

7.30161 

7.30273 

7.30349 

7.30382 

7.30396 

7.30412 

7.30416 

0* 

3 

21.8375 

22.2268 

22.2715 

22.3011 

22.3144 

22.3182 

22.3247 

22.3272 



43.2818 

44.0482 

44.1355 

44.1929 

44.2178 

44.2269 

44.2390 

44.2459 


n 

57.6280 

61.0874 

61.4944 

61.7596 

61.8704 

61.9180 

61.9679 

61.9892 

(b) Potential energy method (upper bound) 

Setting 

angle, 

9 

Mode 

number 

■ 

■ 

■ 

Bj 

n 

■ 

n 

ra 


1 

3.643190 

3.634715 

3.634639 

3.634639 

3.634639 

3.634639 

3.634639 

3.634639 


2 

7.328086 

7.304394 

7.304332 

7.304331 

7.304331 

7.304331 

7.304331 

7.304331 

0* 

3 


23.00084 

22.33244 

22.33243 

22.33158 

22.33158 

22.33158 

22.33158 

4 



45.60155 

44.25278 

44.25276 

44.25100 

44.25099 

44.25099 

44.25099 





5 



66.57146 


62.05877 

62.01729 

62.01712 

62.01699 





oc .uoJ*fo 














































TASLf ill. - EFFECT OF ROTATION UPON THt FREQUENCY RATIOS, p/Ap OF 
PRE TWIST ED, PRECONED BLADES NEGLECTING CORIOLIS EFFECTS 

y « Q, K ■ 0 


TttlCkness 

o/b 

ratio 

Y» 

Degree 

a 

•1 

Mode 

number 

*>C " N* 

Ap c » 30* 

e Pc - 45* 

B Pc - 60* 

•pc * 90* 




1 

3.9400 

3.7697 

3.5576 

3.3111 

3.0450 



0.6 

2 

7.0567 

7.0088 

6.9427 

6.8760 

6.8087 




3 

22.450 

22.356 

22.225 

22.095 

21.964 




1 

5.0*65 

4. *4986 

3.6751 

2.5996 

Unstable 

u .50 

0 

1.0 

2 

7.1263 

6.9371 

6.6669 

6.3849 

6.0899 




3 

23.652 

23.291 

22.789 

22.276 

21.752 




1 

7.7133 

6.0017 

1.9294 

Unstable 

Unstable 



3.0 

2 

11.055 

8.9112 

4.5417 

Unstable 

Unstable 




3 

33.692 

31.578 

28.103 

24.124 

19.345 




1 

3.9516 

3.7925 

3.5637 

3.3189 

3.0548 



0.6 

2 

6.9684 

6.9390 

6.8/09 

6.8022 

6.7328 




3 

22.746 

22.653 

22.526 

22.398 

22.269 




a 

6. 0072 

4.4931 

3.6738 

2.6047 

0.24471 

0.50 

30 

1.0 


7.0755 

6.8788 

6.6005 

6.3097 

6.0049 





23.926 

23.571 

23.078 

22.574 

22.060 




■ 

7.6748 

5.9420 

1.6216 

Unstable 

Unstable 



3.0 

1 

11.064 

8.9271 

4.6029 

Unstable 

Unstable 




B 

34.034 

31.742 

26.305 

24.360 

19.691 

! 



■ 

3.9537 

3.7946 

3.5658 

3.3211 

3.0572 



0.6 


13.189 

13.155 

13.109 

13.062 

13.014 




fl 

24.032 

23.948 

i 

23.633 

23.718 

23.603 




■ 

5.0121 

! 

4.4970 

3.6769 

2.6079 

0.27314 

0.26 

30 

1.0 

I 

13.336 

13.210 

13.031 

12.845 

12.653 




B 

25.099 

24.774 

24.324 

23.861 

23.406 




n 

11.020 

8.6796 

4.502G 

Unstable 

Unstable 



3.0 


14.180 

13.247 

11.822 

10.117 

7.8232 




fl 

34.721 

32.492 

29.180 

25.463 

21.188 





4.5376 

4.1767 

3.6255 

2.9727 

2.1283 





20.556 

20.342 

20.047 

19.748 

19.443 



1 

B 

59.171 

58.976 

58.707 

58.437 

58.165 




n 

5.0131 

4.4979 

3.6776 

2.6089 

0.2816 

0.06 

30 

1.0 

B 

21.049 

20.723 

20.269 

19.804 

19.329 




B 

39.688 

59.387 

58.971 

58.551 

58.127 





11.039 

8.8937 

4.5152 

Unstable 

Unstable 




H 

29.759 

27,720 

25.619 

21.009 

16.574 



i 

H 

69.373 

67.634 

64.336 

60.976 

57.270 















































TABLE IV. - EFFECT OF ROTATION UPON THE FREQUENCY RATIOS, p/ij, OF 
PRECONEO BLAOES INCLUDING CORIOLIS EFFECTS. 


_£L Hwl * *Pc ■ 1& * »Pc ■ *Pc ■ 45 ' PPc - 60' Up c - 90* 


3.9015 3.6346 


2.9719 2. 6551* 


7.3043 7.4956 7. 6607 7.8064 


22.443 22.332 


4.7452 3.8423 2.6358 1.8415 UnsUble 


7.5342 8.1221 8.6401 9.0132 

23.625 23.193 22.599 21.999 


6.5731 3.6629 .55818 2.0659 4.8405 
12.972 14.599 15.694 16.*/o 16.153 
33.533 30.433 26.237 21.871 17.040 


1 3.9039 3.6354 3.2950 2.9718 2.6L66 

2 7.0739 7.2389 7.4311 7.5963 7.7416 

3 22.739 22.628 22.47b 22.324 22.172 


4.7323 3.6276 2.8230 1.8346 

7.4866 8.0744 8.5890 8.9573 


23.898 23.467 22.877 22.263 21.665 


6.5469 3.6354 


.47607 2.0902 4.9093 


12.9698 14.588 15.669 16.114 

33.656 30.547 26.363 22.049 


3.7599 3.5016 3.2331 2 9510 
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TABLE V. - PERCENTAGE VARIATION OF FRE r ' ;NCY RATIOS* DUE TO CORIOLIS EFFECTS 
f < Oi I > 0) L/b * 5, 8p^ * 45 


Thickness 

ratio 

d/b 



Y * 0* 

Y - 30* 

I mode 

11 mode 

III mooe 

1 mode 

11 mode 

111 mode 

-7.379 

-22.837 

-71.070 

7.964 

29.597 

245.553 

-0.207 

-0.834 

-6.640 

-7.540 

-23.159 

-70.642 


-0.222 

-8.710 

-6.861 

-0.250 

-2.200 

-5.851 

-0.004 

-0.064 

-0.208 

0.002 

-0.006 

6.254 

-0.848 

-5.517 

-10.850 

0.671 

4.192 

7.816 

0.041 
0.400 | 
1.679 


‘Percentage variation « (FR C - *"K^ C ; * iuu/i k^ c 
FR Nc - Frequency ratio without Coriolis effects 
Fl£ - Frequency ratio including Coriolis effects 


TABLE VI. - EFFECT OF CORIOLIS TERMS ON THE STABILITY 
BOUNDARY FOR PRETWISTED, PRECONED, ROTATING BLADE, 
MOUNTED OFF THE AXIS OF ROTATION 


d/b - 0.2b, L/d « 20, , = 0, y = 45* 


R 


Pn «= 0* 

p Pc 


8 pc = 15* 6 pc « 45* 


|H 

(a) 

(a) 

2.2 

2.385 

2.464 

2.6 

1.919 

1.962 

3.0 

1.627 

1.648 

4.0 

1.241 

1.238 

5.0 

1.060 

1.038 

8.0 

.859 

.823 

12.0 

.764 

.726 

16.0 

.560 

.543 

20.0 

.411 

.386 

24.0 

.320 

.289 

40.0 

.168 

.127 

100.0 

.060 

.014 


.B23 .823 .348 .347 

.726 .728 .196 .195 


B Pc 

- 60* 

(a) 

(b) 

5.439 

5.414 

3.566 

3.548 

2.390 

2.377 

.841 

.827 

.111 

.107 


.543 .543 .130 .130 

.386 .386 .041 .041 


•Coriolis effects disregarded. 
b Cor1ol1$ effects included. 




















TABLE VII. - EFFECT OF PRETNIST ON THE STABILITY BOUNDARY FOR PRECONED 
ROTATING BLADE MOUNTED OFF THE AXIS OF ROTATION 
d/b • 0.25, i m 0* (Coriolis effects included) 


Value of R for 

Bp c * 15* and 

Value of R for 

Sp c * 45* and 

j- - 3 

*i 

o 

CNJ 

H 

aU" 

r ■ 3 

*i 

f- “ 16 
*1 

1.650 

0.407 

1.871 

0.145 

1.649 

.405 

1.872 

.143 

1.647 

.398 

1.874 

.139 

1.645 

.386 

1.877 

.130 

1.641 

.373 

1.882 

.115 

1.636 

.358 

1.888 

.093 


1.631 


.341 


1.896 


.063 














25 


29 


3 17 

flu x 

(b) d/b *0.05, Pp c - 45°. 


Figure 2. - Effect of variation of pretwist upon the frequency 
parameter ratio of 45° preconed blades at various rotational 
speeds. (Coriolis effects included) 



i3 


OF. POOw Vi» 


d/b-a05,r*30° 

(CORIOLIS EFFECTS INCLUDED) 


PPc , deg 


90YV75 V 0 

i i W i \ i i i i — i — i 

4 .8 1.2 1.6 2.0 2.4 2.8 3.2 3.1 

Qfui 

Fiaure 3. - Effect of variation of precone angle upon the frequency param 
eter ratios at various rotational speeds of a small thickness ratio pre* 
twisted blade. 



Figure 4 - Rotating beam clamped off the axis of rotation 
Pp -p.y-OP. 
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cussed. Furthermore, it Is shown that the inclusion of Coriolis effects is 
necessary for blades of moderate to large thickness ratios while these effects 
are not so Important for small thickness ratio blades. Finally, the results show 
the possibility of buckling due to centrifugal softening terms for large values 
of precone and rotation. 
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